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c§ : ABSTRACT 

^ ' It is well established that certain detached eclipsing binary stars exhibit ap- 

' sidal motions whose values are in disagreement with calculated deviations from 

Keplerian motion based on tidal effects and the general theory of relativity. Al- 
though many theoretical scenarios have been demonstrated to bring calculations 
into line with observations, all have seemed unlikely for various reasons. In par- 
ticular, it has been established that the hypothesis of a third star in an orbit 
almost perpendicular to the orbital plane of the close binary system can explain 
the anomalous motion at least in some cases. The stability of triple star systems 
with highly inclined orbits has been in doubt, however. 

We have found conditions that allow the long-term stability of such systems 
so that the third-body hypothesis now seems a likely resolution of the apsidal 
motion problem. We apply our stability criteria to the cases of AS Cam and DI 
Her and recommend observations at the new Keck interferometer which should 
be able to directly observe the third bodies in these systems. 
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Subject headings: binaries: eclipsing — celestial mechanics — stars: individual 
(AS Camelopardalis, DI Herculis) — instabilities 



Introduction 



Discrepancy between observation and theoretical predictions of the apsidal motion of 
certain detached binary stars has remained an outstanding problem for two decades. In the 
cases of AS Cam and DI Her, for example, observed apsidal motion rates are a fraction of 
the theoretical predictions based on stellar structure, tidal, and relativistic effects. 



As was first pointed out by lRudkj0bing I (119591 ) . the effect of relativistic gravity is 
significant in the case of a number of detached binary stars. Although he considere d only 



DI Her in detail, t he number of interest ing cases has grown t o about half a dozen (IKoch 
I977I : iMoffat lEgsi ). including AS Cam jMalonev et al. IllOSoh . 



The discovery of anomalous apsidal motion by iMartynov and Kha. 



iulhn I fll980f ) was 



initially considered a possible challenge to general relativity. IMoffat I (119841 ) invented an 
alternative theory of gravity that harbored differing predictions for the apsidal motion of 
binary stars and yet maintained agreement with other tests of general relativity. The predic- 
tive power of this theory is weakened by the existence of a new adjusta ble parameter for each 



star. M oreover, increasingly severe tests of general relativity such as in lTaylor and Weisberg 



( 119891 ) make such large deviations at AU scales seem unlikely. Several other, less exotic, so- 
lutions have been proposed. In one scenario the rapid circularization of the orbit occurs 
because of dissipation of angular momentum from stellar oscillations or a large amount of 
mass loss. Another reasonable guess is that the close binary system (CBS) orbit is sur- 
rounded by a resisting medium in the form of gas clouds. The required density well exceeds 
observatio nal limits in the case o f DI H er, however. T hese and other alternatives are re- 



viewed in iGuinan and Malonev I (Il985[ ). iMalonev et al. I (Il989r ). and IClaret I (Il997l . Il998[ ). 



It is the purpose of this paper to demonstrate the feasibility of one particular explanation 
that has been put forward in the literature. Our considera tions may also find application in 
other triple star or suspected triple star systems such as in ICoe et al. I (120021 ) where a triple 
star model for the X-ray pulsar AX J0051-733 is proposed to explain puzzling features of 
the spectrum of a candidate optical counterpart. 

It is well known that the hypothesis of third stars in outer orbit s of these close binary 



systems can bring theory into line with observed apsidal periods (IKhaliuUin et al.l Il991 



Khodvkin and Vedenevev 1119971 ). but the stability of such triple star systems has been in 
doubt (IHarrington I Il968l ). We show here that the inclusion of the apsidal motion as an 
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additional perturbation leads to the conclusion that such triple s tar ra otioi is can be stable- 



In par ticular, the triple star models of iKhodykin and Vedeneyev I (119971 ) and lKhaliuUin et al. 



( Il99ll ) which reconcile the cases of AS Cam and DI Her with observations are shown to be 
stable. Furthermore, we show that observations utilizing the new Keck interferometer should 
be able to directly image the putative third bodies in these two systems. 



2. The Lagrange Planetary Equations for the Hierarchical Three Body System 

We have studied numerically the dynamical evolution of a hierarchical triple system 
consisting of a massive CBS and a third star of moderate mass. Figure [T] shows the nota- 
tion used. Th e calculations were done perturbatively using the disturbing function method 



( iKopal Ill978l ). We have assumed that the three stars are pointlike and isolated from other 
stars. We ignore internal dynamical exchanges such as synchronization, angular momentum 
exchange, and orbital precession of the CBS. Classical tidal effects and relativistic effects are 
assumed to be independent and additive. 



The disturbing functions for the CBS and the third-body are adapted from lBrown &: Shook 



( I1933I . p. 14). They are, respectively. 



^ ^^'^ (mi + ms)" \r'J ^ ^ ' 



1711 + 7712 + ma mim2 



R - — ■ R (2) 

(mi + 7712) 

We average R and Rtb over the mean anomalies of the CBS and third-body (hence, twice 
averaged). The first-order terms of the twice-averaged disturbing functions, R2 and R2, are 
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where the masses are included via the ratios q = 7712/ rrii and q' = m' /rrii. We choose our 
units of measurement to be AU, years, and solar masses, so that the Newtonian gravitational 
constant is G = 47r^. The orientation of the third-body orbital plane with respect to the 
close binary orbital plane is described by the direction cosines {Q,S,N) of the unit vector 
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normal to the third-body orbital plane. We refer the direction cosines to the periastron, a 
perpendicular to the periastron, and the direction normal to the close binary orbit. Let e be 
the angle between the two orbital planes, and call the angle measured from the periastron 
to the line of intersection of the orbital planes (p. We then have 



Q — sine sin (5) 
S — — sine cos (6) 
N = cose (7) 

Since we are averaging over the mean anomaly, there is no Lagrange planetary equa- 
tion for M and furthermore the semimajor axis a has no time dependence. The planetary 
equations for the perturbation of the remaining CBS orbital elements by the third-body are 



where 
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The planetary equations for the perturbation of the third-body orbital elements are 



= 



= 
= B 



(14) 
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dt tb 
, d^' 



where 
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In writing equations ( fT6] - [T8l) we have introduced the direction cosines of the nodal hne 
of the third-body referred to the same CBS axes by which Q,S, and are defined above. 
These cosines are 



F = sin(r2' — Q) cos i sin u + cos(i7' — Q) cos u 
G = sin{Q' — Q) cos i cos u — cos{Q' — Q) sin u 
H = sin(fi' — i7) sini 



(20) 
(21) 
(22) 



Finally, we use cosines of the direction perpendicular to the third body nodal line and 
behind the visual plane. These cosines are 



T = [cos(f2' — i7) COS2 cosi' + sinz sini'] sinu; — sin(fi' — i7) cosi' coscu (23) 
U = [cos(f2' — i7) cosi cosi' + sinz sini'] cosa; + sin(r2' — fi) cosi'sino; (24) 
V = cos(r2' — fi) sinzcosz' — cosisinz' (25) 



3. An Alternative Formulation of the Equations of Motion 

It is worth noting that the equations of motion can be written in a particularly elegant 
fashion by noting that they determine nothing more than an instantaneous rotation that can 
be referred to the CBS coordinates. The probl em is then to express du /dt, di/dt, d(j)/dt, and 
de/dt in terms of this angular velocity vector (iGoldstein et al. II2002I . p. 174). Referring to 
this angular velocity as we first resolve it into three components: ^ = "^^j+^i+^z- These 
components are related by rotation(s) to the angular velocity components. For example, 
is just —di/dt referred to the nodal line Referring to figured! we see that a rotation by 
uj about the z-axis brings this angular velocity into the CBS system; thus 
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cos u sin LJ 
- sin u cos u 
1 



(26) 



where X, Y, and Z are the components of ^ in the CBS coordinates. Proceeding in 
this manner, one obtains three equations giving X, Y, and Z in terms of du/dt, di/dt, and 
dVt/dt. These equations are inverted to yield 



— = cot t[X smuj + Y cos uj\ + Z 

(JjL 

di 

— = —X cos a; + y sin a; 
dt 

— = — CSC i[X sinu + Y cos u] 

(JjL 

Similar considerations for a circular third-body orbit result in 



(27) 
(28) 
(29) 



de 

— = X cos + F sin I 
dt 



dt 



— cot e{X sin (p — Y cos ( 



(30) 
(31) 



Comparing equations (12714291) with (fT0HT2|) . we determine the rate of rotation of the 
CBS frame, as measured by CBS coordinates, to be 



(X, Y, Z) = -A[ ivgii^, NS, + 5Q^ - 2 



(32) 



It has been pointed out by iKiseleva et al. I (119981 ) (using different angles to define orbit 
orientations) that this formulation admits two exact integrals and leads to a first-order 
elliptical differential equation for the eccentricity. Our calculations, however, are entirely 
numerical and are based on equations (ISHT^ and (fTlHTSl) . 



4. The Instability Problem of Hierarchical Triple Star Systems 



There are a couple of ways of seeing the problem of instability. One approach utilizes 
the conservation of angular momentum. We assume that i ~ 7r/2, so that the cotz term in 
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equation ffTOj) can be ignored. This leads to a large eccentricity excursion due to the angular 
momentum exchange between the CBS and the third-body. Since the disturbing function 
R^2 depends on neither M' (the mean anomaly) nor uj' (the longitude of the periastron of 
the third-body with respect to the ascending node), there is no secular variation of the 
third-body semimajor axis a' and eccentricity e'. Therefore, the orbit of the third-body 
maintains its shape, and the magnitude of its orbital angular momentum L' is a constant of 
motion. The direction of the third-body orbital angular momentum, however, will change, 
as the following argument shows. The orbital angular momentum of the close binary system 
Lbs is 3-10 times smaller than L' for the systems under investigation. Since the rotational 
angular momentum of the stars is about 2 orders of magnitude less than L', the total angular 
momentum of the system is Ltot = Lbs + L'. Calculations reveal that the CBS angular 
momentum is transferred to the third star. Since the magnitude of L' can not change, this 
angular momentum transfer forces a change in the orientation of the third-body orbit with 
respect to the total angular momentum. Conservation of angular momentum dictates the 
connection between e and e shown in figure [2l As angular momentum is transferred, the 
coordinate values (e, e) slide along one of the integral curves determined from the initial 
values of the eccentricities e and e', mass ratios q and g', relative inclination e and ratio of 
the semimajor axes a /a'. In the presence of a small but nonzero cotz the above argument 
does not take the orbit all the way to e = 0, but the eccentricity may become close enough 
to unity that tidal effects or even collision destroy the system. This turns out to be the case 
for DI Her and AS Cam. 

We can gain additional insight into the dynamics of the eccentricity by looking directly 
at the equations of motion. More specifically, we look at equation (fTOj) neglecting the cot i 
term and further assume that the motion of the nodal line rjri' caused by nutation and 
precession of the orbits is small. Under these assumptions. 



^ A[2-^Q'-N^] (33) 

at tb 

« -(f) (34) 

at tb at tb 

First, we note that for e < 30° we have Q,S << 1. We see then that the righthand 
side of equation (fTOj) is always positive. Thus u and change monotonically and run over 
all quadrants, never coming to rest. Therefore, the perturbations in the eccentricity (equa- 
tion (^) are periodic, and the eccentricity suffers no excursion to a value nearing unity. 
The situation is much different if e > 30°. At such inclinations, equation ffTOl) implies in 
general four values of (p for which duj/dt = 0. Let us examine the motion of (p at these roots. 
Equations fl33|) and flM|) together imply 
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sm esm 2(p[—) 



(35) 



tb 



dt tb 



This tells us that the four points of stationary u (and 0) will be stable only if sin 20 < 0. 
Looking at equation ([9]) we see that this same condition guarantees {de/dt)tb > 0. Thus 
will wander until it reaches a value that results in an eccentricity excursion. Again, we 
emphasize that these semi-quantitative arguments are limited to the extent that we neglect 
the coti term in equation (JTOll . 



The characteristic time scale for the change of eccentricity can be obtained from equa- 
tion ©: 



For DI Her and AS Cam these times are about 700 and 400 yr respectively, as was 
confirmed directly by numerical integration. Thus at first glance the conclusion seems to 
be that only nearly coplanar hierarchical triple systems can be stable for more than a few 
hundred years. If this conclusion were correct, the third-body hypothesis would be eliminated 
as a probable solution to the apsidal motion discrepancy. 



We now explore the consequences of including the effect of stellar structure (namely 
tidal-rotational deformation of the CBS pair) and the relativistic effect as additive perturba- 
tions on the motion of uj. We assume that the structure effect {duj/dt)ci and the relativistic 
effect {duj/dt)rei act in simple superposition with the effect of the third-body {dij/dt)th, so 
that their influence can be represented by simply adding them to {doj /dt)tb- If the combined 
effect of these two additional terms is of the same order or greater than the third-body effect, 
that is, if 




(36) 



5. A Possible Resolution of the Problem of Instability 




(37) 



then the motion of uj will not stop. Thus oo and will change monotonically, resulting 
in periodic perturbations of the orbital elements of the CBS leading to stability, as in the 
case of low inclinations discussed above. 
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We can derive stability criteria on the basis of equation fl37l) . We consider the cases 
in which either the classical deformation effect or the r elativistic effect don ainates using 
well-k nown relationships for (duj/dt)r i from iKopal I (Il978l ) and {duj/dt)rei from lRudkj0bing 
( 119591 ) and iMartynov and KhaliuUin 



If the classical effect dominates in the sense that \{duj/dt)tb\ < {duj/dt)ci, we have 



If the relativistic effect is predominant [\{duj/dt)th\ < {duj/dt)rei], then 



,2\3/2 



1 



a \ 3 q'a 
a/) Mi(l + g)2 \A 



— cos e I < 5, 



^rel 



G 

72 



10" 



(39) 



Table [U displays the application of these criteria to AS Cam and DI Her. We see that, 
according to our hypothesis, both of these supposed triple star systems are predicted to be 
stable. The stability of AS Cam is provided by the classical effect whereas DI Her is stable 
because of the relativistic apsidal motion. 

We have confirmed this behavior by means of numerical integrations of the equations 
of motion. Th e calculations were done perturbatively using the disturbing function method 
( IKopal Ill978l ). We have assumed that the stars are pointlike, and the close encounters are 
ignored. We ignore internal dynamical exchanges such as synchronization, angular momen- 
tum exchange, and orbital precession of the CBS. The apsidal motion in the CBS caused by 
both the classical tidal-rotational deformation of the components and th e relativistic apsidal 
motion is described by disturbing functions as in iKhaliuUin et al.l (jl99ll ). The classical tidal 
effects and relativistic effects are assumed to be independent and additive. The results for 
AS Cam are shown in figure [2l Thus the angular momentum is transferred back and forth 
between the third-body and the CBS. These periodical variations in e and e can be visualized 
as a flapping of the CBS and third-body orbits, almost like a butterfly (figure [3]) • 



6. Discussion 



We begin our di s cussio n by contrasting our stability criter ia with those of iRoy I (119791 ) . 
Szebehely and Zare I (119771 ). an d lEggleton and Kiseleva I (119951 ) for the case of DI Her using 
the same orbital parameters as iKhaliuUin et al.l (Il99ll ). 
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Roy I (119791 ) allows a very close orbit of the third-body, the restriction being only that 
the semimajor axis satisfy a' > 0.3 AU corresponding to a period P' > 18.3 days. This 
seems much too close to the inner binary for stability at any inclination of the third-body 
orbit. 



Although ISzebehely and Zare I (119771 ) deal primarily with the case of coplanar orbits, 
they indicate how their results may be extended to third-body orbits inclined to the inner 
binary orbit. We have applied their stability criteria assuming that the third-body orbit is 
perpendicular to the inner binary orbital plane. The result is that stability criterion requires 
a' > 2.9 AU {P' > 1.5 yr). This seems more reasonable, but in fa ct numerical simulations 
suggest that this is still too close (IKhodykin and Vedeneyev 1119971 ). 



Finally, lEggleton and Kiseleva I (119951 ) predict that the system is stable for any third- 
body orbit with a' > 1.2 AU (P' > 0.39 yr). 

The criteria we propose prove more restrictive. From equation (l39l) we determine that 
a' > 9.5 AU {P' > 9 yr). We have confirmed the stability of the purported three body system 
of DI Her under this re s triction on the third- b ody o rbit. We also wish to emphasi z e that the 



criteria of iRoy I (Il979h . ISzebehely and Z 



J are 



(119771 ) , and lEggleton and Kiseleva I (119951 ) are 
all based on the purely Newtonian gravitational theory of point mass orbits. Our criteria 
depend on structure effects and/or general relativity. In the case of DI Her the effect of 
general relativity dominates. It appears that the stability of the hypothetical three body 
system of DI Her is to be found in the physics of general relativity rather than structure 
effects. Thus general relativity itself provides stability to the three-body model of DI Her, 
which seems necessary to bring its theoretical apsidal moti on into line with obse rvations. 
A similar stabilizing role for general relativity was found by iHolman et al. I (119971 ) in their 
investigations of a planet orbiting one star of a binary system. Our criteria are somewhat 
more general in that structure effects are also considered. 

We close with a discussion of the prospects of making direct observations of these hypo- 
thetical third-body companions of AS Cam and DI Her in light of recent advances in optical 
interferometry and adaptive optics. 

Indirect evidence for a third-body in AS Ca m ( B-I-B9.5, P=3.43days, e=0 . 17, V =8.6) 
has already been found by iKozyreva et al.l (119991 ) and iKozyreva and KhaliuUin I (119991 ) , who 
found, imposed upon the timing of eclipse minima, a cyclic variation with a period of 2.2 
yr. They have inte rpreted this signal as due t o the Roemer-like influence of a third star. 
The calculations of iKhodykin and Vedeneyev I (119971 ) indicate that in order to account for 
the anomalous apsidal motion of AS Cam, the third-body shou ld be of about IMq. Using 
this estimate, binary masses of 3.3 and 2.5M0 (jHilditch 1972), and the 2.2 yr period, the 
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semimajor axis of the orbit would be 3.2 AU. This corresponds to a hght-t ravel time of 27 



minut es. Combining this with the amphtude of 4.18 minutes measured by iKozyreva et al. 



( 119991 ) yields an orbital inclination of 9° with respect to the line of sight, so that nearly the 
full 3.2 AU is visible to the observer. Assuming a distance of 480 pc, the maximum angular 
elongation is 0".007. This is twice the resolution limit of the Keck interfero meter operating 



at 1.5 /xm with its 85 m baseline. The inte rpretation of eclipse timings by iKozyreva et al. 



(Il999l ) and IKozyreva and KhaliuUin I (119991 ) also predicts the times of maximum elongation. 



In the case of DI Herculis, no indirect indication of third light exists. However, infor- 
mation from past theoretical analyses provides hope that present in terferometers should be 
capable of directly observing the putative third-body. The analysis of iKhaliuUin et al.l (Il99ll ) 
suggests that the minimum third-body mass is about 0.8 Mf7^ and that th e minimum period 
is about 7 yr. Assuming binary masses of 5.15 and A.52Mq (jPopper 1119821 ) . we conclude that 
the semimajor axis is at least 8 AU. Since the orbit is expected to be highly inclined, and 
the distance to DI Her is about 500 pc, we expect a maximum angular elongation of 0".02. 

Of course, the ability to resolve these third bodies is of little use unless they are suffi- 
ciently bright. Here infrared observations are a great advantage, since the compact binary 
stars are relatively massive in comparison with the hypothesized third bodies. For example, 
applying the mass-luminosity relationship to AS Cam, we conclude that the third star should 
have a total luminosity less than the system by 4.3 mag. A simple calculation based on the 
Planck distribution and assumed temperatures of 20, 000 and 6000 K for the binary and third 
star, respectively, predicts that, in the H (1645 ± 155 nm) and K (2200 ± 480 nm) bands, 
the third star is dimmer by only about 1 mag. A similar calculation for DI Her predicts that 
in the H and K bands the third star is dimmer than the system by 3.5 mag. Finally, the 
magnitudes of these systems are such that the compact binary stars may serve as natural 
guide stars for adaptive optics in the case of the Keck interferometer. 

Bolstered by indirect evidence in the case of AS Cam and dynamical stability indicated 
by the considerations of this paper, the case for a third-body solution to the long standing 
problem of anomalous apsidal motion is stronger than ever. The final judgment, however, 
may soon be expected from the current generation of interferometers. 
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Fig. 1. — Kinematic variables describing the relative orientation of the orbits of the CBS 
and the third body. 
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eccentricity 

Fig. 2. — Evolution of e and e due to angular momentum exchange between the CBS and 
the third-body. The loop schematically illustrates the results of numerical integration of the 
equations of motion of AS Cam including the influence of the CBS apsidal motion as an 
additional perturbation. 
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j^i /Ltot 

L', /Ltot 7 



Lbs 



Fig. 3. — Motion schematic. As the CBS loses angular momentum, its orbit overturns and 
eccentricity increases to a value close to unity. If the apsidal motion of the CBS is included 
as an additional perturbation, then e and e oscillations are moderated, and the triple system 
orbits move back and forth like a butterfly. 
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Ta.l)le 1: Stal)ility criteria for AS Cam and DI Her. 
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